Doniach Diagram of Disordered Electron Systems 
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The quantum phase diagram of disordered electron systems (DES) with magnetic impurities (Mis) 
is derived. The competition between the indirect exchange interaction |Jrkky| of magnetic impu- 
rities and the Kondo effect is known to give rise to a rich quantum phase diagram, the Doniach 
DiagrampJ. A Kondo screened phase is separated from a spin ordered phase at a quantum critical 
point (QCP), when the local exchange coupling J and the concentration of magnetic moments n 
are varied. We report here numerical results for disordered 2D electron systems which show that 
both the Kondo temperature Tk and |Jrkky| are widely distributed and the QCP is extended to 
a critical region. Calculating the distribution function of their ratio \Jrkky\/Tk for various sepa- 
rations R, we find a sharp cutoff, and from that a critical density n c (J). For lower densities that 
ratio is smaller than one for all sites, which we can identify therefore as the Kondo phase. As the 
disorder amplitude is increased, a phase of coupled magnetic moments (CM) grows at the expense 
of the Kondo phase. In addition, a paramagnetic (PM) phase where magnetic moments remain 
free at all temperatures arises below a critical exchange coupling J c . Besides a low temperature 
Kondo Fermi liquid phase we find a phase with non-Fermi liquid behaviour with anomalous power 
law temperature dependence in the magnetic susceptibility and the specific heat. We confirm that 
these anomalous powers are related to the multifractality parameter ao- The coupled moment phase 
shows a succession of 3 phases: a Griffiths phase with anomalous power laws determined by the 
distribution of Jrkky, a spin glass phase and a phase with long range magnetic order. We also 
report results on a honeycomb lattice, graphene, where we find that the magnetic phase is more 
stable against Kondo screening, but is more easily destroyed by disorder into a PM phase. 
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Phenomena which emerge from the interplay of strong 
correlations and disorder remain a challenge for con- 
densed matter theory. Magnetic impurities (Mis) in met- 
als are known to stir up the electronic Fermi liquid. These 
correlations result in a strong enhancement of the resis- 
tivity below the Kondo temperature Tk- Impurities are 
also known to result in Anderson localisation and lead to 
an exponential increase of resistivity at low electron den- 
sities. The interplay of Kondo effect with Anderson lo- 
calisation has only recently received increased attention. 
However, spin correlations and disorder effects are rele- 
vant for many materials, including doped semiconductors 
like Si:P close to the metal-insulator transition Q, but 
also typical heavy Fermion systems, like materials with 
4f or 5f atoms, notably Ce, Yb, or U0, 0|. Many of these 
materials show a remarkable magnetic quantum phase 
transition which can be understood by the competition 
between indirect exchange interaction, the Rudcrman- 
Kittcl-Kasuya-Yoshida (RKKY) interaction between lo- 
calised magnetic moments 0-0] and their Kondo screen- 
ing. Thereby, one finds a suppression of long range mag- 
netic order when the exchange coupling J is increased 
and Kondo screening succeeds. This results in a typical 
quantum phase diagram with a quantum critical point 
where the T c of the magnetic phase is vanishing, the Do- 
niach diagram [lj]. In any material there is some degree 
of disorder. In doped semiconductors it arises from the 
random positioning of the dopants themselves, in heavy 
Fermion metals it may arise from structural defects or 



impurities. As noted already early |8|, the physics of ran- 
dom systems is fully described by probability distribu- 
tions, not just averages. This must be particularly true 
for systems with random local MIs[|| . In fact, it has been 
noticed that a wide distribution of Tk of Mis in disor- 
dered host metals gives rise to non- Fermi liquid behavior, 
such as the low temperature power-law divergence of the 
magnetic susceptibility H U^M- It was foundQ 
that non-magnetic disorder quenches the Kondo screen- 
ing of Mis due to Anderson-localisation and the forma- 
tion of local pseudogaps at the Fermi energy. This re- 
sults in bimodal distributions of Tk with a low Tk tail 
and a finite concentration of free, paramagnetic moments 
(PMs) . From this distribution, the anomalous power law 
temperature dependence of the magnetic susceptibility 
and specific heat has been derived analytically and re- 
lated directly to multifractality parameter an [13]. 

Accordingly, in order to be able to understand the 
magnetic properties of disordered magnetic systems, one 
needs to know also the distribution of the RKKY in- 
teraction Jrkky between different Mis 0-0]- ^RKKY is 
mediated by the conduction electrons, and aligns the 
spins of the Mis ferromagnetically or antiferromagnet- 
ically, depending on their distance R. This is a long- 
ranged interaction, with a power law decay that is not 
changed by weak disorder 2lH24| . It is known, that its 
amplitude has a wide log-normal distribution in disor- 
dered metals [HI, [i^]. Thus, to obtain the quantum phase 
diagram, the Doniach diagram of random systems with 
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Mis, we need to derive the full distribution functions of 
both Tk and Jrkky- This is the aim of this letter. 

We present the results of numerical calculations, based 
on the kernel polynomial method, (KPM) [24, 26[ 27 1 and 
compare them with previous analytical results, as ob- 
tained from the multifractal distribution of local inten- 
sities and their correlations fl9| . We start from a mi- 
croscopic description of the Mis, the Anderson impurity 
model coupled to a non-interacting disordered electronic 
Hamiltonian with on-site disorder. Then we map it with 
the Schrieffer- Wolff transformation on a model of Kondo 
impurities coupled to the disordered host electrons (l9|. 
For the numerical calculations we employ the single-band 
Anderson tight-binding model on a square lattice of size L 
and lattice spacing a, 
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(id) 



£F ~ Ecdgc) cjCj, (1) 



where t is the hopping energy, Cj (cj ) annihilates (creates) 
an electron at site i, uii is the on-site random disorder en- 
ergy distributed uniformly in the interval [—W/2,W/2]. 

denotes the nearest-neighbor hopping, e F is the 
Fermi energy as measured from the band edge, e e dge = 
— At and periodic boundary conditions are used. We set 
Planck constant H and lattice spacing a to unity in the 
following. 

We obtain Tk from the Nagaoka-Suhl equation [28| : 



1 = 



D ^ tanh[(e - e F )/2T K ] 
de Poo(£j, 

£ — £p 



(2) 



with band width D. The local density of states (LDOS) 
is defined by poo(s) = (0\S(e — H)\0). J is the local 
coupling constant between the localized MI and conduc- 
tion electrons, sp is Fermi energy. The RKKY coupling 
■^RKKYy between two Mis located at positions Rj, Rj , 
is in the zero temperature limit (T = 0) given by 2J-[26( 



jRKKYi 



-J- 



S(S + 1) 
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de / de 

€<6f Js'>Ef 



,F{e,e')i 



(3) 

where F(e,e')ij = Re[/9y (e)pji(e')], and S is the mag- 
nitude of the MI spin. Using the KPM, one can eval- 
uate the matrix elements of the local density of states 
Pij {e) = (i\S{e - i?)|j>0|2j|2i|3Qi with a Polynomial 
expansion of order M. In order to ensure convergence, 
but avoid finite size effects, we found that M should be in- 
creased linearly with system size L[25], for an explanation 
see Ref. [3li. Eq.© can be solved in a clean 2D system 
analytically @, yielding J° D = sm(2k F R)/(k F R) 2 
in the asymptotic limit (kpR 3> 1). m* = l/(2a 2 t) 
is the effective electron mass, k F the Fermi wave vec- 
tor. A direct comparison between the exact analytical 
expression in the clean limit and the numerical result 
using the KPM method yields excellent agreement. In 




FIG. 1. (Color online) iV( Jrkky) (a) for various distances R 
at W = 2t, (b) for various disorder W at R = 5a (N = 3000, 
L — 100a, M = 1000). Black dashed lines: fit to log- normal 
distribution, (c) P(Tk) for various disorder amplitudes W 
and J/D = 0.25, (d) P(T K ) for W = 5t and various j = J/D 
(N=30 000, L=40a, M=200). ef = It in a)-d). 



accordance with Refs.[2J and [32j, we find that the ge- 
ometrical average amplitude of RKKY interaction de- 
cays exponentially for length scales exceeding localisa- 
tion length £, e ^2 lnJ RKKY> ^ e~ R ^. In Figs.[TJa,b we 
show the distribution function of the absolute value of 
Jrkky for various disorder amplitudes W and distances 
R. Black dashed lines are fitting curves based on the 
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log-normal distribution, P{x) 
where N = 3000 is the total number of disorder configu- 
rations and x = In | Jrkky | ■ The fitting gives for R = 5a 
and W = 2t,At the parameter xq = 5, 6 and the width 
a = 5.3 + .85W/t increasing with the disorder strength 
W. We find that it does not depend strongly on the dis- 



tance R, consistent with analytical results[21[. Next, we 
present in Fig.[T]c the distribution of Tk obtained from 
the numerical solution of Eq. ^ for L = 40, j = 0.25. 
For every sample only one single site is taken to avoid a 
distortion of the distribution due to intersite correlations. 
Therefore, we had to use a huge number of N = 30 000 
different random disorder configurations for each disor- 
der amplitude W to get sufficient statistics. As found 
previously [l3, 14, 33|, the distribution has a strongly bi- 
modal shape where the low Tk- peak becomes more dis- 
tinctive with larger disorder amplitude W. In Fig.Q]d 
we compare these results with an analytical expression 
for the low T^-tail of P(T K ) in 2D DES, obtained in 
Ref.[l9| from the multifractal distribution of intensities: 
In 2D DES there are weak wave function correlations 
which are logarithmic with an amplitude of order 1/g. 
For weak disorder, j> 1, we can rewrite this correlation 
as an effective power law with power r\iD = 2/7rg, where 
g = e f t. The correlation energy in 2D is of the order 
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FIG. 2. (Color online) Critical exchange coupling J c /D (L — 
100) (a) as function of W (M = 300,£f = 2i). Red dashed 
line: Eq. ©. (b) as function of Fermi energy (M = 200). Red 
line: density of states (DOS). 



FIG. 3. (Color online) The probability of the ratio between 
the absolute value of RKKY interaction | Jrkky | and Kondo 
temperature Tk for various distances R as indicated in the 
inset (N = 3000, L = 100a, e F = t,M = 300). 



of the elastic scattering rate, E c 2D ~ 1/t. Thus, we 
find P(T K ) for T K < Max{A c = Z?/£ 2 , A = D/L 2 }^, 
where £ is the localisation length and D = 8t, 



P{T K ) = (1 - p FM ) 



(mn{S,L})-*»S, (4) 



where pfm = ^fm(0)/tim = (Min{£, L})~ 2 "2d is the 
ratio of free PMs, which has to be substracted from the 
weight of P(Tfc)- We see that Eq. (H|) has a power law 
tail with power = 1 — j. 

For T K >T K > Max{A ? = D/£ 2 , A = D/L 2 } we find 
that the power of the tail changes as [l9| 



P(Tk) 
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where c\ = 7.51. Here, £ is the 2D localization length 
in absence of a magnetic field. For parabolic dispersion, 
valid for < cf — e c dgc -C D, one finds £ = <?exp(7rp). 
The scattering rate 1/t is related to disorder amplitude 
W as 1/r = t:W 2 /6D. In Fig.[T]d we compare with 
Eq. (O, using the analytical expression for r/2D, fitting 
E c « .13t and the prefactor. Thus, there is no depen- 
dence on j and the power depends only on t\id- Only 
for the lowest Tk C A, we find a deviation from Eq. 
with a power law that increases with lowering j as 1 — j 
in good agreement with Eq. (fj|. 

Next, we consider the quantum phase transition be- 
tween the freeparamagnetic moment phase (PM) and a 
Kondo phase [ljj. We calculate the critical exchange cou- 
pling J c above which there is no more than one free PM 
in the sample volume L d . From the multifractality of the 
cigenfunction intensities we have found that J c is related 
to the power t\id of the power law correlations in the 2D 
DES as J c = ^/i]2dD. Thus, we find that J c increases in 
2D linearly with disorder strength W asfl9j|. 



Jr = 



D 

3ej 



-W. 



(0) 



In Fig.[5]a, we present the numerical results together with 
Eq. ([5]) as function of disorder strength W and find good 



agreement, with expected deviations at large disorder 
g < 1, where the 1/g expansion breaks down. According 
to Eq. (|6]) J c increases as Ef is moved towards the band 
edge. We have evaluated J c for all Bp and present the 
result along with density of states in Fig.[2]b. We find 
that J c is increasing towards the band edge as 1/i/eF in 
agreement with Eq. ^ . Far outside of e e dge of the clean 
system, it increases as expected with a gap in the DOS, 
as J c /D = l/ln|e edge - e F \- 

Doniach Phase Diagram. In clean systems the crit- 
ical density n c = l/R^ above which the Mis are cou- 
pled with each other is obtained from the condition 
T K . This yields in 2D with 
J 2 s -2^ B 2 and T K = ce F exp(-D/J), 



that U£ KKY (-R C )| 



J 



RKKYlfeFfl»l 



8n 2 k 2 F R 2 

c fa 1.14 that n c = 167r 2 c^ exp(— -j). In disordered 
systems, Tk at a certain site competes with the RKKY 
coupling with another MI at a distance R. Thus, we need 
to consider the distribution function of the ratio of these 
two energy scales x jk = \ Jrkky{R)\/Tk for a given dis- 
ordered sample with density of Mis n = , where R is 
the average distance between the Mis. The distribution 
otxjK for W — it and J/D = 0.2 is shown for several R 
in Fig. [3] We see that it drops sharply at large values of 
xjk , allowing us to define a critical density n c (J) = 1/ R 2 
below which the Kondo effect dominates in the competi- 
tion with RKKY interaction at all sites, see Fig.rj]a. This 
Kondo phase is split at finite temperature into a Kondo 
Fermi-liquid (FL) phase, where all Mis are screened, and 
for T > T*(J) a Kondo Non-Fermi-liquid (NFL) phase 
where some Mis remain unscreened and contribute to the 
magnetic susceptibility with an anomalous temperature 
dependence, given by[19(, 



X(T) 



!2d ( T 
i (tV- 1 ,-, 



for 
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Here, rj = 2(ao — d) depends for d > 2 on the universal 
multifractality parameter ao- In d = 2 it is a function 
of W, V2d{W) = 2/ng(W). The temperature T*(n), see 
Fig-HJc (blue line), is given by the position of the low Tk 
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(a) square lattice 



(b) honeycomb lattice 
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(c) Doniach Diagram 




FIG. 4. (Color online) Critical MI density n c as function of 
J/D for (a) square lattice (ef = t), (b) graphene (£f = 3t, 
Dirac point). L, M, N as in Fig.[3] For n > n c Mis are 
coupled (CM), (c) Schematic Doniach diagram: temperature 
T divided by J 2 versus J/D. Red dot, vertical dashed line: 
critical point J c (n) separating Griffiths phase from Kondo 
phase. Blue line: T*(n) separating Kondo FL phase from 
Kondo NFL phase. For n < and J < J c a paramagnetic 
phase (PM) appears. 



peak in P(Tr-), see Fig.Qjc. 

When the MI density n exceeds n c , magnetic clus- 
ters start to form at some sites and the Mis may be 
coupled by Jrkky(-R)- When R is larger than locali- 
sation length £ the coupling Jrkky(-R) is exponentially 
small and there is a paramagnetic phase (PM) below 
= l/£(g) 2 , where Mis remain free up to exponentially 
small temperatures. For n > n c there is a succession 
of phases, starting with the Griffiths phase where clus- 
ters are formed locally and anomalous power laws are 
observed when the clusters are broken up successively 
as temperature is raised. From the log-normal distribu- 
tion N(Jrkky) one obtains for the magnetic suscepti- 
bility, x{T)T = n FM {T) = J Q T dJ RKKY N(J RKKY ) ~ 
n M exp[~ln 2 (T/J RKKY )/{2a{W)) 2 }, where the width 
a(W) depends on disorder strength W. Accord- 
ingly, the excess specific heat is C(T) = T^§f- ~ 
aip[-ln 2 {T/J RKKY )/(2a(W) 2 }. At higher concentra- 
tions n > nsa a spin-glass phase appears, where the 
magnetic susceptibility shows a peak at spin glass tern- 
Above a 



perature Tsg as studied in Refs. [34 
critical density np a phase with long range order may 
form (35rl38| . In graphene the pseudogap at the Dirac 
point quenches the Kondo effect below J c = D/2, inde- 
pendently on disorder W . Thus, there is a larger param- 
eter space where the Mis are coupled (CM), see Fig. 3] 
(b). However, short range disorder localises the electrons, 



cutting off the RKKY-interaction and for n < there 
is a PM phase. Thus, we conclude that in graphene the 
magnetic phase is more stable against Kondo screening 
but is more easily destroyed by disorder. 
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